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Oh: 

^ : Abstract 

^"T ■ Recently, in an interesting work 'arXiv: 110 6.39721 a solution of the equa- 

[■^^ . tions of motion of massive gravity was discussed, and it was shown that one of 

^D \ the fluctuations on that solution is superluminal. It was also stated that this 

C^ ' rules out massive gravity. Here we find that the solution itself is rather un- 

physical. For this we show that there is another mode on the same background 
which grows and overcomes the background in an arbitrarily short period of 
time, that can be excited by a negligible cost in energy. This solution is trig- 
gered by the parameter governing the superluminality. Furthermore, we also 
show that the solution, if viewed as a perfect fluid, has no rest frame, or that 
C^ \ the Lorentz transformation that is needed to boost to the rest frame is superlu- 

minal itself. The stress-tensor of this fluid has complex eigenvalues, and could 
not be obtained from any physically sensible matter. Moreover, for the same 
setup we find another background solution, fluctuations of which are all stable 
and subluminal. Based on these results, we conclude that the superluminality 



o 



X 



found in larXiv: 1106.3972 1 is an artifact of using an inappropriate background, 
nevertheless, this solution represents an instructive example for understanding 
massive gravity. For instance, on this background the Boulware-Deser ghost 
is absent, even though this may naively appear not to be the case. 



1 Introduction and summary 

Massive gravity (classical theory) has had a turbulent past and present. To briefly 
account for works immediately relevant to the present paper: Fierz and Pauli (FP) 
constructed a ghost-less and tachyon-free linear theory \T]. Van Dam and Veltman, 
and Zakharov (vDVZ), have independently shown [2J that the FP theory has discon- 
tinuity in the zero mass limit, and argued that this excludes massive gravity. Soon 
after, Vainshtein showed that the vDVZ discontinuity is an artifact of the perturba- 
tive expansion that breaks precociously, and argued that upon inclusion of nonlinear 
terms there should be nonperturbative continuity to the massless theory, at least for 
the physical systems of observational relevance, thus evading the vDVZ conclusion 
[3]. However, subsequently Boulware and Deser (BD) [4J showed that in a broad 
class of nonlinear extensions of the FP theory one is not able to retain the needed 
five degrees of freedom of a massive graviton; instead, the sixth mode becomes prop- 
agating on certain backgrounds. This mode typically has negative energies, and is 
referred as the BD ghost. 

More modern developments were triggered by the DGP model [5], for which it 
was argued by Deffayet et.al. [6] that the Vainshtein recovery does take place for 
sources of observational interest. This was followed by a covariant effective field 
theory formulation of massive gravity by Arkani-Hamed, Georgi, and Schwartz [7], 
who also proposed a program to construct a theory that would avoid the sixth mode 
(the BD ghost), starting from the analysis of the decoupling limit [11 El El HO] , where 
things are easier to handle. 

A positive progress toward this goal was made only recently: in Ref. |Tl] it was 
shown that in the decoupling limit the BD ghost can be avoided order-by-order 
to all orders. The absence of the BD ghost in the decoupling limit is a necessary 
consistency condition, but also turned out to be a powerful requirement leading to 
resummation of an infinite number of terms of the effective theory, resulting in a 
covariant Lagrangian with just a few terms [12]. The obtained Lagrangian reads: 

C = ^V^ {R + m\C^l{iq + a,6!l{iq + a.C^ilC))) . (1) 

The tensor /C is defined as follows 

K = ^'i- V^^°9„r9,0V, (2) 

where the square root above denotes a matrix element of the root of the matrix; 
Tjab = diag(— 1, 1, 1, 1), and (p°'{x), a = 0, 1, 2, 3 are four spurious Stiickelberg scalar 
fields introduced as a redundancy to provide for manifestly covariant description of 
massive gravity (for earlier works introducing these scalars, see, [13j.) Finally, the 
mass and potential terms in ([1]) read as follows: 

4ete = ['^']-m\ (3) 

4tl(^) = [JCr-S[)C][)C']+2[}C'], (4) 

/:(^2(/C) = [/C]^-6[/C2][/C]2 + 8[/C=^][/C]+3[/CY-6[/C^], (5) 



where we use the notations [/C] = (Tr /C(;), [Kf = (Tr/C(;)2, while [JC^] = Tr {JC^JC^J. 
The terms Cj'Jj. give total derivatives upon substitution /C(^ — )■ d^d^i:, as indicated 
in their notation. 

The above Lagrangian ([T]) has three free parameters (one of them being the 
graviton mass m), and for some values of these parameters the theory has been shown 
to be free of the BD ghost away from the decoupling limit up to (and including) 
the quartic order in nonlinearities [12]. Remarkably, Hassan and Rosen [14J have 
managed recently to show that it is free of the BD ghost away from the decoupling 
limit, to all order^. We note that the absence of the BD ghost guarantees the 
absence of the sixth mode. This however, does not prohibit one or more of the 
physical 5 polarizations to flip the sign of their kinetic terms on certain backgrounds 
and become ghosts. Such backgrounds should be considered unstable in the theory, 
but such cases should be distinguished from the ones with the sixth mode. For some 
work on cosmology and spherically symmetric solutions in massive gravity see, e.g., 
[T7] - [23], and Ref. [25] for a theory review. 

Recently, in a brief work Gruzinov |26j has found a certain solution of the theory 
([1]), and showed that there is a fluctuation about this solution which is superluminal. 
Based on this observation, it was concluded that massive gravity is ruled out. Below 
we examine this conclusion more carefully. In Section 2 we show that there exist 
a growing solution which overcomes the background arbitrarily quickly, and can be 
excited with virtually no cost in energy. Moreover, this solution is not related to the 
BD ghost, as the latter is absent in this theory. In Section 3 we show that the solution 
of [26j , if interpreted as a perfect fluid, has a stress-tensor with complex eigenvalues. 
Hence, the rest frame for this fluid can only be achieved via superluminal boosts. As 
such, this configuration could not, as an exact solution, be obtained from any known 
physically meaningful form of matter. In the appendix we discuss the decoupling 
limit of the linearized fluctuations. 

The present work does not claim to exclude all possible superluminalities in 
massive gravity. Indeed, some of the terms obtained in the decoupling limit of 
massive gravity resemble the Galileon theories [2Ij, which were shown to exhibit 
superluminalities. It is therefore reasonable to expect that when massive gravity 
reduces to a Galileon theory in the decoupling limit, the fluctuations of the helicity-0 
mode around a spherically symmetric solution can also be superluminal. However, 
for more generic values of the parameters of the theory ([T|), in particular when 
as + 4a4 7^ 0, the decoupling limit of the theory cannot be written explicitly in a 
Galileon form, and it is possible that the decoupling limit does not capture the entire 
physics of the system. In that case a more careful treatment is required, and it is 
yet unclear whether the superluminalities around non-trivial backgrounds survive. 
We plan to report on this issue in a future work. 



^ Our results [12] , as well as the results of [M] , are in conflict with the claim of Ref. [H] of the 
existence of the BD ghost in the quartic order. This controversy is addressed in Ref. [16], where 
it is shown that [151 missed a constraint. See also discussions in Section 2. 



2 Superluminality and growing solutions 

The system of equations of the theory ([1]) reads as follows: 

G^,((?)+m2X^,((7,0) = O, (6) 

m'V^X,,{g,<P) = Q. (7) 

Here, X^^, is a tensor obtained by variation of the mass terms in ([T]), and is given 
explicitly in Section 3. Ref. |26j considers a classical scalar field configuration 



/,2^^^3^ 
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{t,x + et,y,z), (8) 



with an arbitrary constant e, and studies a flat space fluctuation of 0^ in the x 
direction, showing that this fluctuation is superluminal. Based on this, the work 
states that massive gravity is ruled out. 

We start by presenting the results of Ref. [26] in more detail. We do this for a 
small value of e ^ 1, which is enough for our purposes. For this consider the field 
configuration (jHj). It produces some stress-tensor X^^, which also depends on the 
metric g^^; since X^^ is multiplied by m^ in eq. (jH]), one assumes that the back- 
reaction of ni?X^y on the metric is negligible. In this approximation, the remaining 
equation is just an empty space Einstein equation, which certainly has a solution 
g^y = Tj^y. Hence, to summarize the solution of [26] : 

g^lv = Vfiu + h^y, hf,y ~ 0{em^xx-) , (9) 

where for convenience we have included the correction h, where x' denotes some 
components of x^ (these corrections can straightforwardly be calculated for small e, 
and take the form, /zqi oc ern^t^, hi2 oc ern^ty). Since m can be arbitrarily small, 
one can neglect h in ([9]), at least in some region of space and time, and consider 
fluctuations on the approximate background g^u ~ rj^y. 

To present the results of [26] more explicitly, we write down the quadratic La- 
grangian for the fluctuations C,"" of the four components of the 0" field, while freezing 
all the other fields in the theory (the approximation in which this is justified will be 
discussed later, see below and the appendix): 



LO' ^O, 



-"^^°^*^^^ii^- f^") 



The Lagrangian for ("" follows from an expansion of ([T]) on the background ([HI [9]), 
and in the quadratic approximation for the fluctuations C° takes the fornu: 



A = -7^i(C) - ^F.^Fo" - '-Fo,id^C) + 0{e') , (11) 



^Note that C", in spite of its appearance, does not transform as a vector under diffeomorphisms, 
instead, it transforms as a four-coordinate. 



where /3 = —(03 + 1/6), and the total derivative terms have been ignored. Further- 
more, Ref. |2S] focused on the C,2 component propagating only in the x direction 
(this does not excite other fluctuations in (TTTl)). For C,2it)X) the above Lagrangian 
reduces to 

A. = \CI - ^Cl - ^GC'^ + 0{e^) , (12) 

where an over-dot and prime denote t and x derivatives respectively. The dispersion 
relation that follows to leading order in e, w ~ p(l — 3/3e/2), is superluminal, since 
for any nonzero /3, the value of e can always be chosen to give superluminality |26] . 

If the Lagrangian (fT2!) is taken in isolation of all the other fields and interactions, 
as done in [26] , then, the derived superluminality can be removed by a simple change 
of coordinates to x and t where, x = x + t(3/3e/2) and t = t (which is just a galilean 
transformation with velocity equal to — 3/3e/2). Likewise, from an innocent field 
theory of a scalar ip coupled to a source J with the Lagrangian, —{d^ip)'^ + ipJ, 
one can get the Lagrangian of the type flT2l) by the above change of coordinates. 
Hence, if flT2|) were the entire Lagrangian one could quantize fluctuations in the 
{x, t} coordinate system where no superluminality would appear. 

The actual question, however, is what and how these fluctuations couple to other 
fluctuations and external sources, and what those other fluctuations do. In the full 
theory the field C, does mix with the tensor and scalar modes at the linearized level, 
and has also nonlinear interactions. Also, there are 0{e^) terms neglected in (TT2|i ^l. 
We study in turn all the fluctuations omitted in (fT2l) . 



Let us first focus on other components of C" which were not considered in 
Dropping the last term in (ITTl) . we write the equations of motion in the Lorentz 
gauge d^C'' = 0: 

□(Co-gCi) + g5o5iCo-g5oCi = o, 

niCi + qCo) + qdodiCi-qdlCo = 0, 
na + qi-did,Co + 2d,doCb - ^o^ftCi) = , (13) 

where b = 2,3, and q = 3(3e/2. These empty-space equations have many growing 
solutions. The one we focus on is 

Ci - -qml e + mlt{x - x^) + O (e^) , (3 ^ -ml t{y - y^) + O (e^) , (14) 

where rriQ, xq, yo, are arbitrary integration constants, and other components of ( are 
set to zero. In the leading order in e, the e^ pieces in the above expressions should 
be ignored. Note that for flT^ . Fqi = qmlt + ml{x — Xq), and F02 = —ml{y — y^). 



■^ We thank Mehrdad Mirbabayi who pointed out to us that the Lagrangian (fTT|) , if considered 
in isolation, has a gauge symmetry, since the last term in it can simply be removed to the next, e^ 
order, by a field redefinition, Co ~^ Co ~ fCij Ci ~^ Ci ~ fCo- This symmetry, is not present in the 
e^ order, and also in the full massive theory, due to the coupling of Ca to the tensor mode. 



At each point in space one could choose corresponding Xo,yo, such that Fqi = qmlt 
and Fq2 = point-by-point in the whole space. 

There are some important comments to be made about the solution. 

(a) The solution flT^ grows in time t and overcomes the background ([8]) for 
t > t*, where t* = (mMpi/mg). However, the mere existence of this solution cannot 
be interpreted as an instability of the background (|8]). The reason is that the solution 
dH]) has a nonzero energy density proportional to Fq^ + F^^, and in order to excite 
this field configuration within a finite volume in space one would need some energy; 
if such an energy is supplied, it is then not surprising that (fT^ does overcome the 
background after some period of time. 

Nevertheless, there is an aspect of the solution ( !T^ that suggests that the back- 
ground (E]) is unphysical. To see this let us introduce a length scale L^ = qt^,. Let 
us now consider a small imaginary box of volume L^ centered around the point 
X = xq, y = yo, z = 0. Most importantly, we take L ^ L^,. Furthermore, imagine 
that we supplied enough energy in the box, and the appropriate boundary condi- 
tions at its sides, so that inside the box the solution flT4|) is excited, while outside of 
the box the background is still given by ([S]). Let us now estimate how much energy 
density we need to supply for this to be the case. The average energy density in the 
box will consist of three terms, D = D1 + D2 + D3, where Di ~ q^m^t^ , D2 ~ qtm^L, 
and D^ ~ ttIqL^. Now, for any time moment t > L/q, the Di term dominates. What 
is important, however, is that L/q ^ L*/g = t*. Hence, at the time moment t ~ t*, 
when the solution (fT4l) in the box begins to dominate over the background (IE]), the 
energy density that is required to excite it is of the order Di ~ g^; however, the 
latter happens to be zero in our approximation since we are ignoring terms of order 
e^ ~ g^ in the action. Therefore, we conclude that at the expense of the energy 
density that is zero in our approximation, we can excite the solution flT^ in a vol- 
ume of size L^ ^ Ll, and this solution overcomes the background ([8]) after t r^ t^, 
time. Since Xq and yo iii dSD are arbitrary, we can now consider the whole space 
populated by non-overlapping boxes in each of which an appropriate value of the 
parameters xo,yo are chosen, and as a result, the growth described above develops 
in each of these boxes. Then, it is logical to interpret t* as a characteristic time for 
the growing solution to dominate in the entire space. Most importantly, this time 
scale is independent of q, and can be made arbitrarily small by adjusting the inte- 
gration constant ruQ. Moreover, the amount of the energy density needed to excite 
this solution, Di ~ q'^rriQt'^, although negligible in our approximation, in any event 
is much smaller than the characteristic scale of the stress-tensor for the background, 
em'^Mpi, as long as t < {mMpi/{eml))] the latter is always the case, for t ^ t*. 

(b) Although, the above growth can be arbitrarily fast, it is triggered by the same 
parameter g ~ e that sets the superluminality of the (2 niode. This phenomenon 
disappears in the limit g — )■ 0, even though for g = the solution (fT^ still grows and 
it may appear that the above-given argument would still hold in the g — )• limit. In 
this limit the solution (fT^ reduces to the terms mlt{x — Xq) and mlt{y — yo), but 



no t^ term is remaining. Then, within a box of size L^ (which is now necessarily 
larger than L\ that tends to zero) one would be able to excite the solution at the 
cost of a non-negligible energy density D ~ 1)3 ~ m^L?'. This energy density cannot 
be ignored as there are no neglected q^ terms in the action any more. Clearly, such 
a solution does not indicate any problem of the background, it simply reflects the 
fact that the background has changed at the expense of the supplied finite energy 
density Q 

(c) The Lagrangian (fTTI) . as was pointed out above, is not the total Lagrangian 
of our theory, even at the linearized level and even in the leading order in e. Then 
the question arises of whether the mixing with other fields plays any role. Since this 
is a bit technical, we address this question in the appendix, where we show that only 
in the ?ti — )■ limit the dynamics of the C,a modes in (TTT]) can be decoupled from the 
tensor and scalar modes, assuming that one ignores nonlinear interactions as well. 

In spite of the above described issue with the background ([9]) it is important to 
emphasize that this problem is not related in any way to the Boulware-Deser mode, 
i.e., to a potential sixth degree of freedom in a broad class of massive gravities, which 
is absent in the present model ([T]). To see that the sixth mode is not propagating 
on the background considered in [26], we look at the full action for the tensor fields. 
In unitary gauge 0° = x°, the solution considered in [26] amounts to taking a 
background solution for the metric which is Minkowski, but not in Cartesian form. 
Specifically the background metric in unitary gauge is ds"^ = —dt^ + {dx — edt)"^ + 
dy'^ + dz^. Expanding to quadratic order in perturbations around this solution and 
to first order in e, the mass term which is now expressed entirely in terms of the 
tensor field since C" = is 

2 

Tfl — — — — — — 

- -J- m. -h^ + 6cie/ii,/i° - (6ci + lyhoih) , (15) 

with ci = —as — 3/2. We can use this form in order to count the physical degrees of 
freedom. In order for the BD ghost to be absent, one has to have the Hamiltonian 
constraint [1]. In the FP linearized theory the Hamiltonian constraint is enforced 
by hoo being a Lagrange multiplier, while h^j is algebraically determined by an 
equation that is independent of /iqo- In the Lagrangian (ITB]) . however, /iqo mixes 
with hoi, and this may seem to forbid the presence of a constraint. However, this 
is not so, there still exists a linear combination of the fields that is a Lagrange 
multiplier in the approximation used. A convenient way to see this is to calculate 
the determinant of the 4x4 Hessian matrix for the Lagrangian 'H^''' = ,,^ f, — . If 
the determinant is zero, then there are constraints. It is straightforward to calculate 



^This is similar to the case of a free massless scalar field which has a solution, <j) ^ t; in any 
infinitesimal region of space one needs to supply a nonzero energy density to excite this solution. 
Likewise, to excite the solutions, (j) ~ tx, ty, which also exist in this theory, one would need finite 
energy density in any finite volume. 



that the determinant of the Hessian that follows from (TT5l) is of order e^, i.e., it 
is zero in our approximation, while the rank of the Hessian is 3. Hence, there is 
one constraint in the system. Moreover, conservation of this constraint leads to a 
secondary constraint, as shown for these theories exactly in [T3]; due to these one is 
able to eliminate the BD ghost. 

One may also consider the counting of degrees of freedom in the non-unitary 
gauge considered in [2S] • In such a gauge the tensor mode propagates two degrees of 
freedom, while there should be only three degrees of freedom for the four Stiickelberg 
fields. The latter requirement at first sight seems unlikely since there is no gauge 
invariance for ( in the full theory, and moreover, (o enters with a time derivative, even 
in the simplest case of ci = 0. Based on this one may be tempted to conclude that the 
Lagrangian for the ( field (after gauge fixing h^^) propagates 4 degrees of freedom. 
However, by more careful inspection one can show that there are constraints that 
render only 3 degrees of freedom in the ( sector in general (detailed discussions of 
how this works in the full nonlinear theory are given in [T6]). 

In conclusion, the solution iQ seems problematic, despite being ghost-free. Is 
there another problem- free solution for the very same configuration of the Stiickelberg 
fields ([H])? The answer is positive. It is straightforward to find another solution to 
the system of eqs. (Q and ([7]), for given ([8]): 

9,. = V,u + e{6X + Sy.) + 0{e') . (16) 

or exact to all orders ds"^ = —dt^ + {dx + e dt)"^ + dy"^ + dz'^. The above solution 
differs from (jH]), by e, i.e., by the same parameter that sets superluminality found 
in [26]. Furthermore, it is easy to notice that the solution (IHfTB!) is nothing but 
the Minkowski solution, g^i, = t]^^, cjf = (t,x,y,x), transformed by the coordinate 
change x'^ — )■ x^ + et(5f . Therefore, the fluctuations above the solution (IHfTBl) are 
just ordinary fluctuations of the Fierz-Pauli theory, 

2 

- —r'^v"^ [KuKp - h^^Kp) , (17) 

o 

which are known to be subluminal and stable. 



3 Superluminality of the source 

The solution considered in [26] is not an exact solution of massive gravity. As we 
have explained, it is at best a solution valid locally in a space-time region whose 
size/time scale is set by L ~ l/{y/em). Alternatively we can allow it to be an exact 
solution by adding an external source T'^ which is chosen so that 

m^X^, = T-*. (18) 

In principle we could imagine this external source being set up by a configuration 
of matter, a fluid, or a set of scalar or gauge fields. However it is easy to see that 



the 'fluid' needed would itself be composed of superluminal matter. To see this, 
imagine T^^* were described by a perfect fluid. Let us assume that the fluid has a 
rest frame. If this is the case we can perform a Lorentz transformation so that in the 
vicinity of one point the fluid has zero velocity. At that point T^^^ = 0. Since the 
background metric is flat, goi = and so Tg'^t j = T^^to — 0- This in turn implies that 
the energy density T^oext is one of the eigenvalues of the stress energy tensor T^^ext- 
The stress-energy tensor is expressed in terms of the tensor /C^ in the combination 

X^, = ICg,, - ]C,, + (1 + Sag) U% - /C/C^, + ^ ([/C]^ - [/C^]) g^^ (19) 

+a (kI - ICKl + i/C,, {[ICr - [/C^]) - Im' - 3[/C][/C2] + 2[}C'])g,.^ , 



where a 
Xk 



of JC^l, 



as + ia^. So the eigenvalues of T'^j^ext are determined by the eigenvalues 
which in turn are expressed as 



A 



(n) 



1 



1-A 



(n) 
Y 



for n 



,4, 



where Ay are the eigenvalues of Y^i,: 

r^ = g^'^d^rduc^'va 



(20) 



(21) 



It is straightforward to show that the first two eigenvalues of this tensor are complex 
for the background solution, and so T°oext is complex in this frame. Explicitly for 
the background considered the matrix fl^Tl) is 



/ 1 _ e2 _e \ 
e 10 

10 

\ 1/ 



f22) 



and its eigenvalues are easily shown to be 



,(1) ,(2) ,(3) ,(4) 
Ay , Ay y '^Y ^ ^Y 



1 



-eV4 - e2, 1 - -e^ - -eV4 - e^ 1, 1 . (23) 



We can therefore immediately infer that the eigenvalues of /C are also complex, and 
so are the eigenvalues of T^^cxt- 

This implies that there is no rest frame for the fluid, or that the Lorentz trans- 
formation needed to boost to the rest frame is superluminal (and hence a complex 
transformation), since it is not possible to perform a real Lorentz transformation to 
set Tq^^ = 0. As such this configuration could not, as an exact solution, be obtained 
from any known physically sensible form of matter. 

Even in the absence of a source, the same arguments hold. It is clear that for any 
solution of the equations G^jy + m?X^i, = which looks locally like flat space-time 



8 



with the field profile described in ( 18|9|) . it is not possible to boost to a frame in 
which G°* = in the local vicinity of a point. In this sense the solution already at 
the level of the background looks superluminal and rather unphysical. 

We would like to thank L. Berezhiani, G. Chkareuli, S. Dubovsky, D. Pirtskhalava 
and R. A. Rosen for useful discussions and comments, and especially M. Mirbabayi 
for his valuable input. CdR is supported by the Swiss NSF and GG is supported 
by the NSF grant PHY-0758032. GG regrets that discussions on the subject of 
could not be contained within NYU. 



Appendix 

Below we show that the Lagrangian ( ITTI) can be obtained from the full theory in the 
limit 771 — !■ 0. For this we start with the mass terms on the background (IHlEl) 

Cm = ^ {h% -h^ + e,h,X - ^2hoih + 4eV) , (24) 

where ei = e(6/3— 4) and €2 = e(6/3 — 3), and all the indices are contracted by r]^^. We 
express the Lagrangian in terms of the Stiickelberg fields by using the substitution 



hf^u _ S^^ _ df,CdX T]ab 

Mpi mMpi 771"^ M^y 



V -> tF - ^^ - "' 7..2 ^ (25) 



where we defined S^i, = d^(u + d^fj, + e^S'^dXi + S^d^d), and introduced canonical 
normalizations for all fields. 

Then, the total Lagrangian reads as follows: 

Tn 
C = CEH{h) + C^{0 + -S^''{K.-v,uh) 

+ ^(/iiA^ + 5iX)-^(V5 + 5oi/i) + 0(m2), (26) 

o o 

where CEH{h) denotes the linearized Einstein-Hilbert term, while >Cc(C) is the La- 
grangian given in ( TTTl) . We also ignored the terms of order 0{m^) or smaller, 0{em?), 
in (126|) . Note that the tadpole appearing in fl2^ gets canceled in fl26l) by the corre- 
sponding tadpole coming from the EH term taken on the background ([9]). 

We see that the fields C,"" mix to the tensor field. Our goal is to show that this 
mixing disappears in the m — )■ limit. For this we note that we can remove the 
third term in Cc_{C,) by a linear field redefinition (see footnote 3), and then introduce 
the helicity-0 field by the change of variables C^ ^ Cm + d^ir/m. As a result we get 
the following Lagrangian: 

£ = CEH{h) + C,{C) + ^V^''{7r){h,,-v,uh) + 

+ ^ {h^d^'don + did^nh^) - ^ {hoiDn + dod^nh) + C(m, em, m^) , (27) 



where Ci{C) is the Lagrangian (TTTj) less the last term, and V^uij^) = 28^8^71 + 
^{5^^8y8in -\- 5^8^8iTi) + |(5^9,^(9ovr + 5^c?^c?o7r). As we see, there is a mixing between 
the tensor mode and the helicity-0 mode vr. Due to this mixing the helicity-0 gets a 
kinetic term via the shift h^j, — > h^i, + rj^yTt] as a result, the helicity-0 would couple 
to an external source had we introduced it in the theory. However, the field C, does 
not couple with anybody in this limit. This would be so even if we were to introduce 
a stress-tensor of an external matter. The coupling of C, appears only at a nonlinear 
level. 

References 



[1] 

[2] 

[3] 
[4] 

[6; 

[7] 
[8] 

[9 

[lo; 
[11 

[12 

[13 
[14 



M. Fierz and W. Pauh, Proc. Roy. Soc. Lond. A 173, 211 (1939). 

H. van Dam and M. J. G. Veltman, Nucl. Phys. B 22, 397 (1970); 

V. I. Zakharov, JETP Lett. 12 (1970) 312 [Pisma Zh. Eksp. Teor. Fiz. 12 (1970) 

447]. 

A. I. Vainshtein, Phys. Lett. B 39, 393 (1972); 

D. G. Boulware and S. Deser, Phys. Rev. D 6, 3368 (1972). 

G. R. Dvah, G. Gabadadze and M. Porrati, Phys. Lett. B 485, 208 (2000) 
|arXiv:hep-tE70005016i . 

C. Deffayet, G. R. Dvali, G. Gabadadze and A. I. Vainshtein, Phys. Rev. D 65, 
044026 (2002) |arXiv:hep-th /0106001|. 

N. Arkani-Hamed, H. Georgi and M. D. Schwartz, Annals Phys. 305, 96 (2003). 

M. A. Luty, M. Porrati and R. Rattazzi, JHEP 0309, 029 (2003) 
[arXiv:hep-th/0303116 |. 

P. Creminelli, A. Nicolis, M. Papucci and E. Trincherini, JHEP 0509, 003 

(2005). 

C. Deffayet and J. W. Rombouts, Phys. Rev. D 72, 044003 (2005) 
tarXiv:gr-qc/0505134] . 

C. de Rham and G. Gabadadze, Phys. Rev. D 82, 044020 (2010) 
jarXiv:1007.0443l [hep-th]]. 

C. de Rham, G. Gabadadze and A. J. ToUey, Phys. Rev. Lett. 106, 231101 
(2011); [arXivilOll. 12321 [hep-th]]. 



W. Siegel, Phys. Rev. D 49, 4144 (1994) |arXiv:hep-th/9312117] . 
S. F. Hassan and R. A. Rosen. I^XTv: 1106. 33441 [hep-th]. 



10 



[15] A. H. Chamseddine and V. Mukhanov, larXiv: 1106.58681 [hep-th]. 

[16] C. de Rham, G. Gabadadze, and A. J. Tolley, "Ghost free Massive Gravity in 
the Stiickelberg language" larXiv: 1107.38201 

[17] C de Rham, G. Gabadadze, L. Heisenberg and D. Pirtskhalava, Phys. Rev. D 
83, 103516 (2011) |arXiv:1010.1780l [hep-th]]. 

[18] K. Koyama, G. Niz and G. Tasinato. arXiv: 1103.47081 [hep-th]. 

[19] T. M. Nieuwenhuizen. l arXiv:1103.5912l [gr-qc]. 

[20] S. F. Hassan and R. A. Rosen. larXTv: 1103.60551 [hep-th]. 

[21] K. Koyama, G. Niz and G. Tasinato. arXiv:1104.2143l [hep-th]. 

[22] G. Chkareuh and D. Pirtskhalava. larXiv: 11 05 l783l [hep-th]. 

[23] A. Gruzinov and M. Mirbabayi, larXiv: 1106. 25511 [hep-th]. 

[24] C. de Rham, L. Heisenberg, [arXiv:1106.3312 l [hep-th]]. 

[25] K. Hinterbichler,larXiYillQL3735| [hep-th]. 

[26] A. Gruzinov. larXiv: 1106.3972 1 [hep-th], 

[27] A. Nicohs, R. Rattazzi and E. Trincherini, Phys. Rev. D 79, 064036 (2009) 
[arXiv:0811.2197 [hep-th]]. 



11 



